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In this paper we consider a numerical method for proving a class of analytical 
inequalities via minimax rational approximations. All numerical calculations in 
this paper are given by Maple computer program. 



1. Some particular inequalities 

In this section we prove two new inequalities given in Theorem 1.2 and Theorem 
1.10. While proving these theorems we use a method for inequalities of the following 
form f(x) > 0, for the continues function / : [a, b] — ► R. 

1.1. Let us consider some inequalities for the gamma function which is defined by 
the integral: 

oo 

(1) r w~J e " ttz ~ ldt 



which converges for Re(z) > 0. In the paper [10] the following statement is proved. 
Lemma 1.1 For x € [0, 1] the following inequalities are true: 

(2) T(x+ 1) < x 2 - -x+ % 

v ' 4 5 

and 

(3) (x + 2)T(x + 1) > \. 

o 

The previous statement (Lemma 4.1. of the paper [10]) is proved by the approx- 
imative formula for the gamma function Y{x + 1) by the polynomial of the fifth 
order: 

(4) P 5 (x) = -0.1010678a; 5 +0.4245549a: 4 -0.69985882: 3 +0.9512363a; 2 -0.57486462:+l 

which has the numerical bound of the absolute error e = 5 • 10~ 5 for values of 
argument x£ [0, 1] [3] (formula 6.1.35., page 257.). 
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In the Maple computer program we use numapprox package [6] for obtaining the 
minimax rational approximation R(x) = P m {x)/Q n {x) of the continuous function 
f(x) over segment [a, b] (m is the degree of the polynomial P m (%) an d n is the 
degree of the polynomial Q n {x))- Let e(x) = f(x) — R(x) be the error function 
of an approximation over segment [a, b}. Numerical computation of R(x) given by 
Maple command: 

(5) R := minimax ( f(x), x = a. .6, [m, n], 'err') ; 

The result of the previous command is the minimax rational approximation R(x) 
and an estimate for the value of the minimax norm of s(x) as the number err 
(computation is realized without the weight function). With the Maple minimax 
command a realization of the Remez algorithm is given [29] , [30] . If it is not possible 
to determine minimax approximation in Maple program there appears a message 
that it is necessary to increase decimal degrees. 

Let us assume that for the function f(x) the minimax rational approximation R(x) 
is determined. Then, in the Maple the same estimate for the minimax norm of the 
error function e(x) is given by command: 

(6) err:— infnorm ( e(x), x — a..b ) ; 

The result of the previous command is number err := max l£ [ a ^] \f{x) — R(x)\. 
Practically for the bound of the absolute error function \e(x)\ we use e — err. 
Let us remark that the bound of the absolute error e is a numerical bound in the 
sense [16] (approximation errors on page 4.), see also [17]. 

Let us notice, as it is emphasized by the Remark 4.2 of the paper [10], that for 
the proof of Lemma 1.1. it is possible to use other polynomial approximations 
(of lower degree) of the functions T(x + 1/2) and T(x + 1) for values x £ [0, 1]. 
That idea is implemented in the next statement for the Kurepa's function which is 
defined by the integral: 



(7) K (z) = J , 



t t z ~ 1 , 



which converges for Re(z) > [4]. It is possible to make an analytical continuation 
of the Kurepa's function K{z) to the meromorphic function with simple poles at 
z = — 1 and z = —n (n > 3). Practically for computation values of the Kurepa's 
function we use the following formula: 

= Ei(i) + z7r + (-i)*r(i + z)r(-z,-i) 



which is cited in [9]. In the previous formula Ei(z) and T(z, a) are the exponential 
integral and the incomplete gamma function, respectively. Let us numerical prove 
the following statement: 
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Theorem 1.2 For x € [0, 1] the following inequality is true: 

(9) K(x) < K'{0) x, 

where K (0) = 1.432 205 735 . . . is the best possible constant. 

Proof. Let us define the function f(x) — K (0) x — K(x) for x <G [0,1]. Let us 
prove f(x) > for x e [0, 1]. Let us consider the continuous function: 



(10) g(x) 



a : x = 0, 



for constant a = — ^ — . Let us notice that the constant: 
2 

K"(0) ,. K'(0)-K'(x) ,. if'(0) a: - AYs) ,. /(x) 

(11) a = ^ = lim — — — = lim — — — — = lim ^-V-. 

2 2:^0+ 2x x^0+ x Z x^0+ X 

is determined in sense that g(x) is a continuous function over segment [0, 1]. The 
numerical value of that constant is: 

°° 2 

(12) a = -- I e"*-^-^ dt = 0.963 321 189 .. . (> 0) . 

2 J t 1 

o 

Using Maple we determine the minimax rational approximation for the function 
g(x) by the polynomial of the first order: 

(13) Pi (x) = -0.531 115 454 x + 0.921 004 887 

which has the bound of the absolute error E\ — 0.04232 for values x e [0, 1]. The 
following is true: 

(14) g(x) - (Pi (a;) - £i) > and Pi (a;) - £i > 0, 

for values x <G [0, 1]. Hence, for x e [0, 1] it is true that g(x) > 0, and f(x) > as 
well. Q.E.D. 

Remark 1.3 Numerical values of constants K (0) and K (0) are determined by 
Maple program. The numerical value of K (0) was first determined by D. Slavic 
in [5]. 

Corollary 1.4 A. Petojevic in [12] used an auxiliary result K{x) < 9/5 x, for values 
xG[0, 1], from [10] (Lemma 4.3.) 7 for proving new inequalities for the Kurepa's 
function. Based on the previous theorem, all appropriate inequalities from [12] can 
be improved with a simple change of fraction 9/5 with constant K (0). 
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1.2. D. S. Mitrinovic considered in [1] the lower bound of the arc sin function, 
which belongs to R. E. Shafcr. Namely, the following statement is true. 

Theorem 1.5 For < x < 1 the following inequalities are true: 

3x 6(VTT^-V1~^) 

(15) < ; < arc sin x . 

2 + v/i - x 2 ~ 4 + vTT~? + - 

A. M. Fink proved the following statement in paper [7]. 
Theorem 1.6 For < x < 1 the following inequalities are true: 

3X 7TT 

(16) < arcsmx < 



2 + Vl-x 2 ~ '' " 2 + Vl - x 2 ' 

B.J. Malesevic proved the following statement in [8]. 
Theorem 1.7 For < x < 1 i/ie following inequalities are true: 

TT 

(17) <arcsina:< ^ — 1 < 



2 + VI -x 2 2 + ,/l - x 2 2 + VT^I? 

tt - 2 

Remark 1.8 The upper bound of the arc sin - function: 

TT 

(18) 4>(x) = - *~ 2 ' 



-x 



tt-2 

is determined in paper [8] by X-method Mitrinovic- Vasic [1]. 

L. Zhu proved the following statement in [11]. 

Theorem 1.9 For x e [0, 1] the following inequalities are true: 



3x ^ 6(VTTx-VT^x) ^ 

< , , < arc sin x 



2 + Vi - x 2 A + VTT^ + VT 

(19) i 

tt(V2+-)(VT+x~-VT^) nx 

< , ; < ; 

4 + Vl + z + Vl-z 2 + VI - x 2 

In this paper we give an improved statement of L. Zhu. Let us numerical prove the 
following statement: 
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Theorem 1.10 For x <E [0, 1] the following inequalities are true: 

3x 6(\/TT~F - IF) 

< , , < arc sin a; 



2 + VT 



4 + VI + x + \J\-x 
tt(2 - v/2) 



(20) 



< 



7T -2v/2 



(VTT^-VT^) 



v/2(4 - tt) 

tt - 2\/2 



7r(v / 2+x)(Vl + a;- VI -a:) 
< 1 < 



4 + VT+^+VT T7 ^ " 2 + Vl - a; 2 
Proof. Inequality: 

n(V2 + h(VTT^- VT^) ~VtS (VI + ^ - VI - a) 



(21) 



> 



4 + VT+^ + VT^ ■ ^(4-;) +VTT ^ + vr 



for a; G [0,1], is directly verifiable by algebraic manipulations. Let us define the 
following function: 



(22) 



^-^ (VT+^-VT^) 
r(x) = — = arc sin a;, 

^i^ + VTT^ + VT^ 

7T - 2V2 



for x e [0, 1]. Let us prove /(ar) > for x e [0, 1], ic. f(smt) > for t e [0, -] 
Let us define the function: 



(23) 



g(t) = { 



t = 0, 



/(Si,,,) :«€(0,f) 



/3 



f = 



2' 



where a and /? are constants determined with limits: 



(24) 



and 



a = lim 



/(sin*) 



(4 + V2)7r - 12V2 
(24 - 12V2)tt 2 



> 



(25) = lim A™** = (16V2-16) + (8-4V2K-V^ > 0. 

t-7r/2- /3 fZ[_/\ (2V2-2)7T a 
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The previously determined function g(t) is continuous over [0, — ] . Using Maple 

we determine the minimax rational approximation for the function g(t) by the 
polynomial of the first order: 

(26) P x (t) = 0.000 410 754 1 + 0.000 543 606 

which has the bound of the absolute error E\ = 1.408 • 10~ 5 for values t£ [0, ^] . It 
is true: 

(27) g(t) - (Pi(t) - ei) > and P^t) - El > 0, 

for values ie[0,— ]. Hence, for ie[0, — ] is true that g(t) > and therefore 
/(sini) > for te [0, |] . Finally /(a;) > for a; G [0, 1]. Q.E.D. 

Remark 1.11 The paper [15] considers the upper bound of the arc sin - function: 

?r(2 - V2) , 

(28) (^(x) = ^ - 2 ^ 



■{s/l + x-s/l-x) 



v/2(4 - tt) 

7T-2V2 

via \-method Mitrinovic-Vasic [1]. 



2. A numerical method for proving inequalities 

In this section we expose a numerical method for proving inequalities in following 
form: 

(29) f(x) > 

for the continuous function / : [a, b] — > R. Let us assume that x — a is the root 
of the order n and x = b is the root of the order m of the function /(x) (if x = a 
is not the root then we determine that n = 0, ie. if £ = 6 is not the root then we 
determine that m = 0). The method is based on the first assumption that there 
exist finite and nonzero limits: 

(30) n lim — and = lim ■ /( ''' ) 



z^a+ ( x - a)™ (6 - x) m ^^6- (x - a) n (b - x) m ' 

If for the function f(x) (over extended domain of [a, b} ) at the point x = a there is 
an approximation of the function by Taylor polynomial of n-th order and at point 
x = b there is an approximation of the function by Taylor polynomial of m-th order, 
then: 

f(")(a) f (m Hb) 
(31) a= * W and = (-l) m 777 l \ ■ 

v ; n!(6-a) m v y m! (6 - a) n 
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Let us define the function: 



a 



x = a 



(32) 



9{x) = gl b {x) = < 



x e (a, b) 



v 



(x - a) n (b-x) 

(3 



m 



x = b 



which is continuous over segment [a, b]. For proving inequality (29) we use the 
equivalence: 



which is true for all values x G [a, b]. Thus if a < or [3 < the inequality (29) 
is not true. Hence, we consider only the case a > and (3 > 0. Let us notice 
that if the function f(x) has only roots at some end-points of the segment [a,b], 
then (33) becomes f(x) > iff g(x) > for x G [a, b]. The second assumption 
of the method is that there is the minimax (polynomial) rational approximation 
R(x) = P m (x)/Q n (x), of the function g(x) over [a, b], which has the bound of the 
absolute error e > such that: 



for x G [a, b}. Then g(x) > 0, for x G [a, b}. Finally, on the basis (33), we can 
conclude that f(x) > 0, for x G [a, b]. 

Let us emphasize that the minimax (polynomial) rational approximation of the 
function g{x) over [a, b], can be computed by Remcz algorithm (via Maple minimax 
function [6]). For applying Remez algorithm to the function g(x) it is sufficient that 
the function is continuous. If g{x) is differentiable function than the second Remez 
algorithm is applicable [2]. According to the previous consideration, the problem 
of proving inequality (29), in some cases, becomes a problem of existence of the 
minimax (polynomial) rational approximation R{x) for g = g ab (x) function with 
the bound of the absolute error e > such that (34) is true. Let us notice that 
the problem of verification of inequality (34) reduces to boolean combination of the 
polynomial inequalities. 

Let us consider practical usages of the previously described method. For the func- 
tion of one variable the previous method can be applied to the inequality the infinite 
interval using the appropriate substitute variable, which transforms inequality to 
the new one over the finite interval. Next, if some limits in (30) arc infinite, then in 
some cases, the initial inequality can be transformed, by the means of appropriate 
substitute variable, to the case when the both limits in (30) are finite and nonzero. 

The advantage of described method is that for the function f(x) we don't have to 
use some regularities concerning derivatives. Besides, the present method enables 
us to obtain computer-assisted proofs of appropriate inequalities, which have been 
published in Journals which consider these topics. With this method we have 
obtained numerical proofs of the appropriate inequalities from the following articles 



(33) 




(34) 



R(x) -e > 0, 



[18], [19], [20] [21], [22], [23], [24], [25], [26], [27], [28]. 
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Finally, let us emphasize that the mentioned method can be extended and applied 
to inequalities for multivariate functions by the means of appropriate multivariate 
minimax rational approximations. 
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